NOTICE 


THIS DOCUMENT HAS BEEN REPRODUCED FROM 
MICROFICHE. ALTHOUGH IT IS RECOGNIZED THAT 
CERTAIN PORTIONS ARE ILLEGIBLE, IT IS BEING RELEASED 
IN THE INTEREST OF MAKING AVAILABLE AS MUCH 
INFORMATION AS POSSIBLE 







DOE/NASA/1040-17 
NASA TM-61 569 


PRELIMINARY RESULTS FROM A 
FOUR-WORKING SPACE, DOUBLE 
ACTING PISTON, STIRLING 
ENGINE CONTROLS MODEL 


(N&S&*TH«81569) PBBLIHIMABT RESOLTS FBOH A 
PODR-ROBKIMG SPACE, DOOBLE-ACT ING PISTON, 

STIBLING ENGINE CONTROLS MODEL (NASA) 17 p 
HC A02/HP A01 CSCL 20D 

G3/3U 


Carl J. Daniele and Carl F. Lorenzo 
National Aeronautics and Space Administration 
Lewis Research Center 


Work performed for 

U.S. DEPARTMENT OF ENERGY 
Conservation and Solar Energy 
Office of Transportation Programs 


Prepared for 

Fifteenth Intersociety Energy 
Conversion Engineering Conference 
Seattle, Washington, August 18-22, 1980 



/ 


N30-29o24 
Unci j 3 

2 8*» t> 4 


) 

) 



/ 


NOTICE 


This report was prepared to document work sponsored by the 
United States Government. Neither the United States nor its agent, 
the United States Department of Energy, nor any Federal employees, 
nor any of their contractors, subcontractors or their employees, 
makes any warranty, express or implied, or assumes any legal lia- 
bility or responsibility for the accuracy, completeness, or use- 
fulness of any information, apparatus, product or process dis- 
closed, or represents that its use would not infringe privately 
owned rights. 


DOE/NASA/1 040-1 7 
NASA TM-81 569 


PRELIMINARY RESULTS rROM A 
FOUR-WORKING SPACE, DOUBLE- 
ACTING PISTON, STIRLING 
ENGINE CONTROLS MODEL 


Carl J. Daniele and Carl F. Lorenzo 
National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135 


Work performed for 

U.S. DEPARTMENT OF ENERGY 

Conservation and Solar Energy 

Office of Transportation Programs 

Washington, D.C. 20545 

Under Interagency Agreement EF-77-A-31-1040 


Fifteenth Intersociety Energy 
Conversion Engineering Conference 
Seattle, Washington, August 18-22, 1980 


PKEUMlNAkY RESULTS FROM A FOUR-MOMKINb SPACE. UUUttLE-ACTlNG 

PISTON. STIRLING ENGINE CONTROLS MOOEL* 

by Carl j. Daniele ana Carl F. Lorenio 
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Lends Research Center 
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SUhl«iARY 

A four-Morking space, double-acting piston. Stirling engine simulation 
is being oevelopeo for controls stuoies. Tne development methoo is to con- 
struct tMO simulations, one for detailed fluid oehavior. ana a second model 
with simple fluid benaviour but containing tne four working space aspects 
ana engine inertias, validate these models separately, then upgrade the four 
working space moael by incorporating tne detailed fluid behaviour model for 
all four working spaces. 

The single working space (SwS) model contains the detailed fluid dynam- 
ics. It has seven control volumes in which continuity, energy, ano pressure 
loss effects are simulated. Comparison of the SwS model with experimental 
data Shows reasonable agreement in net power versus speed characteristics 
for various mean pressure levels in the working space. 

The four working space (FwSJ mooel was built to Observe tne oenaviour 
of the Whole engine. Tne drive dynamics ano vehicle inertia effects are 
simulated. To reuuce calculation time, only three volumes are used in each 
working space and the gas temperatures are fixed ino energy equation). 
Comparison of the FwS mooel predicted power with experimental data shows 
reasonable agreement. Since all four working spaces are simulated, tne 
unique capabilities of the mooel are exercised to loox at working fluiu su»/- 
ply transients, short circuit transients, ano piston ring leakage effects. 
The FWS model has been upgraded by using the detailed SwS mooel fur each of 
the four working spaces. Currently the detailed FWS model is being reworked 
to reduce the amount of calculation time per cycle. 


NUmENCLATUKE 


A area, m2 

Cp specific heat at constant pressure, J/(kg-K) 

Cy specific heat at constant volunte, J/(kg-K) 

F force, N 

f( j function of 

G gear ratio 

h neat transfer coefficient, J/(sec-m2-K) 

I inertia, N-m-sec^ 

J mecnanical equivalent of neat, 1.0 (N-m)/u 


*Work funded by the U.S. Department of Energy unoer Interagency Agreement 
EC-77-A-3 1-1040. 


k resistance! (N>sec)/K 9 

t piston stroke, m 

P pressure, Him 

Q heat flow rate, J/sec 

R gas constant, (N-m)/(kg-K j 

R|, connecting rod length, m 

fluid resistance, (N-sec)/(kg-m^) 
T temperature, K 

t time, sec 

^ torque, N-m 

V volume, m^ 

Vel vehicle velocity, Km/nr 

M mass, kg 

A mass flow rate, kg /sec 

Y ratio of specific heats 

a Change In 

e angular position, rad 

Subscripts: 


a auxiliaries 

c cooler 

d drag 

e engine 

f mechanical friction 

h heater wall 

1 position In working space, inaex 

In Into 

j,k working space 

£ piston ring leakage 

m mesh 

0 oeao volume 

out out of 

p piston 

r rod 

rr rolling resistance 

s stored 

SHT snort circuit 

SHTMAX Short circuit maxiiiium 
SHTwiN Short circuit minimum 
sup supply 

volume gas volume 
w wheel 


INTRODUCTION 

A four-working space, oouole-actlng piston, Stirling engine simulation 
Is being developed at NASA Lewis Kesearcn Center for controls stuoles. The 
development methoa Is to construct two simulations wnich together form a 
complete model. The models are: a single worxlng space mooel, with oe- 

talled fluid oenavlour, and a second model, with simple fluid behaviour out 
containing the four working space aspects ana the engine mechanical Iner- 
tias. The approach Is to tnen combine them Into a single uetaileo controls 
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model. The sIngle-MoricIng space (SWS) inoael was developed to determine the 
number of control volumes (lumps^ ano equation types neeoeo to aoequately 
represent the thermouynanilcs of the engine. The SWS model has oeen written 
to approximate a detailed performance mouel (Kef. 1) developed by Tew. Tne 
SWS mooel differs from the performance model In that the SWS mooel has fewer 
control volumes* Includes Integrated flow pressure drop effects ano uses 
average values of neat transfer coefficients and flow resistances determined 
from the performance model results. Also the SkS mooel does not contain 
Shuttle or conduction losses since these only Influence efficiency but oo 
not Impact control behaviour. 

A second mode) was constructed to study the behaviour of tne total en- 
gine system. Tne mooel consists of four working spaces (FwS) between four 
pistons. To reduce the calculation tinm on the digital computer, only tnree 
control volumes are used In each working space and gas teinperatures are 
fixed within each volume (no energy equation). Tne drive dynamics are also 
Included and tne FwS mouel can be run oy either forcing tne piston motion or 
as an engine. 

Included In this paper are the results obtained from botn models com- 
pared with some experimental oata. Also, since all four working spaces are 
simulated* the unique capability of the FMS mooel Is exercised to study var- 
ious phen(Miena not easily duplicated by a single working space model. To 
demonstrate this capability, working fluid supply transients, short circuit 
(braking) effects, and piston ring leakage pertomiance are presented. 

Finally, the capabilities of botn models nave been combined Into a sin- 
gle detailed controls simulation whlcn will eventually be used with a de- 
tailed engine control system simulation to explore various controls concepts. 


SlNbLl: WUKKll^b SPACE (SwS) NUIJEL 

A single working space mouel of a Stirling engine was developed to 
determine an adequate representation for tne fluid dynamics and thermo- 
dynamics of the engine. Seven control volumes were selected one each for the 
expansion space ano heater, three for the regenerator to give a reasonable 
temperature distribution, and one each for the cooler ano compression 
spaces. A schematic of tne mooel Is snown In Fig. 1 


SwS hooel Equations 


The equations useo to mooel the fluid dynanilcs and tnermooynamics are 
slnif I If Icatlons to the complex partial olfferentlal flow equations. Mass 
flOM Is assumed to occur oue to pressure differential between gas mooes. 
From Fig. 1, some representative equations are, for mass flow: 


Wi ■ ^ (Pj - P i+j^) 1 • ... 6 


The Change In stored mass In a volume Is: 




1 « 1 , 1 : 

*0 • *'7 


. . 1 
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The energy equation used to calculate the change In tenq)erature In a volume 
Is complex due to the oscillatory nature of the Stirling cycle. The form of 
the equation (Ref. 2 ) Is; 


**s^^1 " •l-lG^I-l ’ ^l) ‘ *l(^^1 



1 - 1.2 . . . 
Wq ■ Wy ■ U.O 


where the neat flows to ana from the gas are modeled as: 


l3j 


Q - nA(aT) (4) 

The primed varlaoles are Interface volume temperatures and are determined by 
uirnlnd differencing (Ref. 1). The method assumes that the gas temperatures 
have the same temperature profile as tne regenerator matrix temperature. A 
representative equation Is: 



The regenerators are modeled as thermal lags: 


h.A. 

’’m, - (’ 

•^m *m 


'1 



1 - 3 , 4,5 


(0) 


Pressure Is calculated using the Ideal gas law based on volume tempera- 
tures. Variable volumes are calculated using the drive geometry to calcu- 
late piston position as a function of crank angle. These equations are 
described later. 

Equations of the forms of (2), (3), and (6) for the seven volumes re- 
sult In a 17th order model. The Integration technique useo Is a backward 
difference method which utilizes a multivariable Newton-Khapson iteration 
(Ref. 3). This technique was chosen because of the large difference In time 
constants which results when flow resistance Is considered. 

Comparison of the SUS model with the Tew performance model shows agree- 
ment of pressures within 2 percent and temperatures within 8 percent over a 
cycle at steady state. This Is considered adequate for controls analysis 
where performance detail Is not so Important. 
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FUUK UUKKlNIi space (FUS) hUDEL 

The second model generated was a tour worKiny space mooel. Since the 
purpose of the mouel was to Investigate tne four working space aspects of 
the whole engine, the fluid dynamic moael of each working space was simpli- 
fied. Only three control volumes are used In each working space. One vol- 
ume for the heater (and expansion) space, one volume for tne regenerator, 
and one volume for the cooler (ano compression) space. All volumes are as- 
sumed Isotherms l-w1th the expansion ana compression space gas temperatures 
being the average steady state temperatures from the SwS model. The re- 
generator gas and mesh temperatures are assumed to be the same. Mhlle the 
mooel Is quite simple It will oe shown that It gives a fairly good represen- 
tation of the engine. The model also has the orlve dynamics Incluaea as 
well as a simple mooel for vehicle load effects. A schematic of tne model 
Is shown In Fig. 2. 


Fhb Mooel Equations 


Tne fluid dynamic equations used In tne FimS mooel are slinpllfleu from 
tne SWS model. From Fig. 2, the pistons are nuinbereo In tne order In wnicn 
they reach top stroke as Is Inolcateo by the arrows In the piston heads. 
Double subscripted variables Indicate; first subscript, a voluirm position 
within a working space; secono subscript, the working space. Some represen- 
tative equations are: 


'*’l,j * ^^1#j ~ *’l+l.J^/*^1.J»1'^l.j 


Vi.j ■ "'-j ■ 


J • 1.4 

1 « XfZ 
J « lt4 
1*193 




> 




o.u J 


(b) 


and 





J - i.A 
1 ■ 1,3 


Variable volumes are calculated as a function of piston position: 


(b) 


''1,1 * * (2 ■ ^l)^p 

Piston position Is a function of crank angle: 

X, - Q) cos e, ♦ -(l/ Sin^e, - - (|) 


llO) 


1,2 ... 4 (11) 
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Torque Is calculetea as a function of differential forces on the pistons 
which are summeo tnrougn the drive geometry: 





Also includea in the FwS model are simplified vehicle load effects ano en- 
gine power losses due to mechanical friction ano auxiliaries. Figure 3 
shows a schematic of now these losses are calculated. Torques from the four 
pistons are summed to form indicated torque. Torque due to engine friction 
(^) is subtracted to form brake torque. This is then available for aux- 
iliaries (^) and also for vehicle load effects such as rolling resistance 
and drag The vehicle inertia is brought into the torque 

equation through the gear ratio to give an effective engine inertia. The 
summation of torques is integrated to give engine speed ano integrated again 
to give crank angle, which in turn is useu to generate piston position 
(Eq. The torque equations for tne losses are linear and nonlinear 

functions of e with tne torque equation: 




CUFiPAkl^ON UF FwS AND SW^ MODELS WlTh EXHERIMENTmL DATA 

Both the SmS ana FwS moueis were run at various engine speeds from loDu 
to 4000 rpm and three different mean pressure levels ranging from 6 to 
15 MPa. Power was calculated for each cycle until steady state was 
reached. Power calculations were maue by integrating the area of tne 
pressure-volume curves in the expansion and compression spaces; ano then 
adding the results at tne end of the cycle. The power output preuictea oy 
the SlkS ntodel was multiplied oy four to give gross engine power. Tne FwS 
model sums up tne predicted power (indicated) from all four working spaces. 
From these values, mechanical ana auxiliary power losses are subtracted to 
get net power. The results are snown in Fig. 4. both models snow reason- 
able agreement over tne wnole map with tne SWS model snowing slightly better 
agreement than the FwS model at the higher mean pressures ano speeds, 
although neither model preuictea tne fall in power at 4UUU rpm ano iS mPa. 


UNiguE CAPAolLlTItS UF Tht FwS hUDEL 

The purpose of the moueis discussed in this report is for controls 
studies. The FwS model is simplifieo but can oe useo to stuu^ overall en- 
gine perfonnance in response to various control schemes ana inputs. Also 
since all four working spaces are simulated, the FWS mouel can oe useu to 
study phenomena not easily duplicated by a single working space mouel. Some 
applications of these will now be snown. 
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Working Fluid Supply Transient 


The niean pressure control system for Stirling engines modulates engine 
power by changing the amount of working fluid In the cycle. In order to 
Increase engine speed, working fluid Is supplied to the engine. From 
Fig. 2, fluid Is supplied through the piston roos; when the pistons are near 
bottom (mln^ stroke. The fluid Is supplied to the compression space through 
a timing slot in tne piston roo. A block diagram of the supply system Is 
shown In Fig. 5. The cneck valves are open as long as tne supply pressure 
Is greater than tne corresponding compression space pressure. The timing 
slot area varies as a function of piston position. Flow Is modeled as a 
linear function of pressure drop: 

*sup^ * Fj* ^%up " ^3,1^ 1 ■ 1,2 ... 4 (14) 

A typical supply transient Is snown In Fig. 6. Tne starting engine 
speed Is 2UU0 rpm. After tne simulation settled out, working fluid Is 
added. This Is shown as a function of crank angle In tne top grapn. Note 
that the amount of working fluid added decreases with each cycle since tne 
working space pressure level continues to rise. Only the supply for tne 
second working space Is snown. ^Supplies for tne otner working spaces are 
similar but phase shifted by <<U* with some overlap. Note tne Increase In 
torque amplitude when injection begins. This is due to the timing effect of 
the Injections. Also, as fluid Is added, tne compression space pressure 
rises. The mean pressure shown is an average value for tne sinuslooal com- 
pression space pressure. Tne pressure rise Is 3 hPa in u.24 seconds. 

Tne torque wave shape snows four cycles for every engine cycle before 
tne supply transient starts. Tnis occurs due to the summation of tne tor- 
ques from all four pistons. When tne working fluid is aoded to the cycle, 
the torque starts to rise. There Is no transient torque drop seen In the 
engine. This Is probably cue to tne absence of tne energy equation for the 
Incoming fluid in the compression space. The addition of cold working fluid 
to the compression space in tne actual engine has a quenching effect on tr« 
ten^erature rousing pressure to drop wnlcn results In a transient crop In 
torque. Note tnat tne torque continues to show tne effect of tne four work- 
ing spaces during the supply transient by the changes in wave shape during 
tne rise and fall of tne torque wave shape. These changes occur wnen tne 
working fluid Is added to tne various working spaces. Trie rather slow speed 
response can be attributed to tne assumed vehicle weight causing tne effec- 
tive Inertia to be quite large. For the figure snown, engine speed In- 
creased 30 rpm In u.24 seconds (8 cycles} which resulted in a venicle speeo 
Change from 31.86 to 32.5 km/nr. 


Piston Ring Leakage 

The next effect studied is piston ring leakage. Tnis allows working 
fluid to pass from one working space to an adjacent one. This analysis 
would be difficult to duplicate on a single working space model since tne 
flow can redistribute through all four working spaces. It Is assumed in 
this analysis tnat all four piston rings nave the same leakage area. Tne 
leakage terms are designated as At in Fig. 2. The leakage flow is 
modeled using an orifice equation: 
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^/W. , - P, J 

r ‘^'1.3 - ' 3 .^ 
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1»2 . . . 4 
J ♦ 1 

j ■ 4* k • 1 
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Tn«_resuUs art snown In Fig, 7. The nominal leakage area chosen Mas 
0.002 caK which corresponds to a gap of about u.uuu2& cm. The resultant 
effective leakage flow resistance at 15 MPa mean pressure and 4U00 rpm is 
about 30 times greater than the flow resistance in the working spaces. This 
results in a maximum flow rate 13u times larger in the working spaces tnan 
through the leakage area. The maximum flow rate ratio crops with speeo to 
30 to 1000 rpm and IS MPa; and with leakage area increase, to 34 at iUOO rpm 
and 0,000 cm leakage area, Tne simulation was run at four different engine 
speeds from looo to 4000 rpm and six different leak.'ge areas from 0 to 
0.01 cm2. The indicated power was taken when steady state was reached ano 
then the mechanical ano auxiliary power losses were subtracted from it as 
was done in Fig. 4. Mote that the A ■ 0 case is the same as plotted on 
Fig. 4. As tne leakage flow area increased, the resultaht net power of tne 
engine dropped. For a leakage area of 0.004 cm^, a power reouction of 12 
percent at 40o0 rpm and almost 34 percent at 1000 rpm is shown. This analy- 
sis shows the value of a four working space mouel to determine power losses 
that can occur cue to faulty piston rings In a Stirling engine. Tne effects 
of a failure of a single ring can also oe easily Investigated with a four 
working space muoel of this type. 


Short Circuit Analysis 

Rapid engine oraxlng of tne Stirling engine Is accomplished by snort 
circuiting of the working spaces. Tne analysis of a snort circuit control 
mooe Is roost effectively done with a four working space engine rouuel. Fig- 
ure b shows a schematic of the snort circuit system. Tne compression space 
of each working space Is connected to two olfferent plenums one containing 
high pressure (PsHTMAXi containing a low pressure (PsHTMIN^ 

Check valves allow flow to the high pressure plenum when the corresponolng 
compression space pressure is greater than P^>HThAX« A second set of 
check valves let flow to tne compression space wnen tne corresponding c(xii- 
presslon space pressure 1$ lower than PshtkIN* working fluid Is 

removed from the compression space when tne pressure Is at a maximum tnereoy 
lowering the pressure; and added to the compression space when the pressure 
is at a minimum thereby Increasing tne pressure. The system thus oecreases 
the pressure swing In the working spaces ano thus lowers tne power output of 
the engine. A valve connects tne two plenums and determines tne amount of 
engine braking that is accomplished, with the valve Shut, the plenum pres- 
sures reach their respective maximum and minimum values, with tne valve 
open, pressure In the nigh pressure plenum falls ano pressure In the low 
pressure plenum rises. Note that the system requires working flulo to be 
shuttled to and from two different plenums anu four working spaces. This 
would very difficult to simulate using only a single working space mooel. 

Figure 9 shows the net steady state power as a function of speeu ano 
valve opening ratio. Note tne large drop In power as the valve is open for 
any constant speed. For exatik)le. If the valve Is naif open, the power crops 
from 34.5 to 1.25 kw at 4UUU rpm. If tne valve Is open sufficiently, the 
power goes negative which Indicates the engine Is absorbing power from the 



loao. The aate tor these curves t«as ootalneo oy forcing piston position as 
a function of time (orlve oynamlcs not active) ana letting the simulation 
settle out. At the nignest valve openings (u.b to 1.0) the pOMer starts to 
Increase as speed drops. 

The torque curves corresponding to the power curves In Fig. g are shown 
In Fig. lu. Note that tne torque also drops significantly with valve open- 
ing at constant speed. At 4UU0 rpm. It drops from 86 to 0 N-m at u.S area 
ratio. As speed decreases, torque continues to orup off for tne nlgner 
valve area settings. Figures 8 and lo Indicate that snort circuiting Is a 
very effective means of supplying engine oraalng torque. 


60 M^lAKY 

Presented In this report are preliminary results from a simplified four 
working space Stirling engine simulation. The mooel was constructed to sim- 
ulate overall engine oehavlour, rather than precise performance, for con- 
trols studies. While the model Is quite simple It exhiolts reasonaole re- 
sults When compared to steady-state engine test results, and gives ootn 
steady-state and transient results which are representative of Stirling en- 
gine oenavlour. Also presented Is a single working space model which has a 
detailed model of the fluid dynamic ana thermodynamic aspects of a single 
working space. This mooel also shows reasonaole results when compareo to 
experimental data. Currently, the FwS mooel has oeen upgraoeo oy Incorpo- 
rating tne oetalleo fluio oenavlour mooel for all four working spaces. Tnis 
oetalleo model will be a reference mooel for Staling engine controls stud- 
ies. Tne simplified models will oe noniially used when they yield valid re- 
sults. But, in some cases, such as snort circuiting, a detailed four work- 
ing space iiiodel is required. Results from tne detailed mooel should compare 
even more favorably with the Stirling engine oenavlour since tne neat trans- 
fer Characteristics of trie oifferent phenomena simulated can oe incorporated 
Into the model. For example, the quenching effect of the cold supply fluid 
on the compression space temperature can oe included. Tnis r.ould result in 
a better prediction of the torque waveshape curing a supply transient. 

Future efforts will include conceptual controls studies such as variable 
stroke control and detailed modeling of the current mean pressure control. 
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